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Disorder is an unavoidable ingredient of real systems. Spatial disorder generates Griffiths phases 
(GPs) which, in analogy to critical points, are characterized by a slow relaxation of the order 
parameter and divergences of quantities such as the susceptibility. However, these singularities 
appear in an extended region of the parameter space and not just at a (critical) point, i.e. there is 
generic scale invariance. Here, we study the effects of temporal disorder, focusing on systems with 
absorbing states. We show that for dimensions d > 2 there are Temporal Griffiths phases (TGPs) 
characterized by generic power-law spatial scaling and generic divergences of the susceptibility. 
TGPs turn out to be a counterpart of GPs, but with space and time playing reversed roles. TGPs 
constitute a unifying concept, shedding light on the non-trivial effects of temporal disorder. 
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Quenched disorder affects the behavior of particle sys- 
tems, altering critical properties and introducing new 
universality classes. Non-magnetic impurities in mag- 
netic systems or defects in type-I superconductors are 
typical examples of this [lj]. Moreover, novel phases, 
with phenomenology unheard-of in pure systems, can be 
induced by spatial disorder. This is the case of Grif- 
fiths phases (GPs) appearing in classical, quantum, and 
non-equilibrium disordered systems [2j-y|. GPs, which 
are of relevance in condensed matter physics as well as 
in other contexts 0, Q , are regions of the phase space 
—actually a sub-region of the disordered phase —charac- 
terized by extended singularities of the thermodynamic 
potentials and, as a consequence, generic divergences of 
magnitudes such as the susceptibility @]. Furthermore, 
GPs are characterized by an anomalously slow (power- 
law or stretched exponential) relaxation to zero of the 
order parameter (and of other time-dependent quanti- 
ties) which contrasts with the fast (exponential) decay 
typical of pure systems. Such an anomalous relaxation 
in the disordered phase occurs owing to the presence of 
rare regions where the disorder is such that the system is 
locally in its ordered phase and, hence, a potential barrier 
has to be overcome for it to relax. The convolution of dif- 
ferent, exponentially rare, sizes with exponentially large 
decaying times gives rise to an overall slowing down of the 
system's dynamics, which typically becomes algebraic in 
time —with continuously varying exponents —and loga- 
rithmic at the critical point (see below for more details) 
d, H| . Divergences in the potentials together with slow 
relaxation are two features strongly reminiscent of criti- 
cality and its concomitant scale-invariance. However, in 
GPs these traits appear not just at a critical point but in 
a broad extended region, providing a robust mechanism 
to justify some cases of scale-invariance in Nature Q. 



or ecology requires parameters to be disordered in time 
rather than in space Q ■ This is the case of magnetic sys- 
tems under a fluctuating external field [l| , or of ecological 
populations under changing environmental conditions 0] . 
In general, temporal fluctuations in the parameters can 
shift critical points Q and affect universal features both 
in equilibrium ficl ] and non-equilibrium systems (llj . In 
a pioneering work, Leigh showed that, in one- variable 
(mean-field) models of stochastic populations, environ- 
mental noise changes the system mean lifetime (time to 
reach the absorbing state) from exponential to a power- 
law in system-size 0, This result inspired us to 
systematically explore the role of temporal disorder in 
spatially extended systems (beyond mean-field) and to 
study if rare temporal regions induce new phases analo- 
gous to spatial-disorder induced GPs. Do temporal Grif- 
fiths phases exist? If so, which properties do they have? 
Do they exhibit any type of generic scale invariance? 

To tackle these questions, we start by analyzing a spe- 
cific model with absorbing states: the contact process 
(CP) in the presence of temporal disorder. In the 
CP, each site of a d-dimensional lattice can be either oc- 
cupied z(x) = 1 (active) or vacant z(x) = 0. At each time 
step, an active site is randomly chosen and, with proba- 
bility b, it converts into active a nearest neighboring site 
(provided it was empty), while with probability 1 — 6 it is 
declared empty. Time, t, is then increased by 1/N(t), 
where N(t) is the total number of active sites. The 
"pure" CP is critical only at some dimension-dependent 
value b c . vure (d) separating an active from an absorbing 
phase (see |13| and the schematic diagram in Fig[T])- This 
phase transition, occurring at b c . pure w 0.767, 0.622, and 
0.5 for d = 1, d = 2 and d = oo, respectively, lies in 
the very robust directed percolation universality class [l3| . 
For the spatially-disordered case, b is replaced by 6(x); in 
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point of this quenched version, b q:C > b C:Pure (sec Fig. [TJ). 

Temporal disorder is implemented by allowing b to be 
a time-dependent random variable, b — > b(t), for all x. In 
the simplest (uncorrelated) case, b takes a random value 
extracted at each Monte Carlo step (i.e. whenever the in- 
teger value of t increases) from some distribution of mean 
bo and width a. Correlated fluctuations can also be im- 
plemented, by allowing b to obey an Ornstein-Uhlenbeck 
dynamics [lil, EH- This temporally-disordered contact 
process (TD-CP) is similar to the Jensen's model in 
. Following the instantaneous value of b, the system 
shifts between the tendencies to be active (b(t) > b CiPUre ) 
or absorbing (b(t) < b CtPure ), provided that the disorder 
distribution is broad enough (see FigfT] for a schematic 
diagram). Owing to fluctuations, any finite system is, 
however, condemned to end up in the absorbing state, 
either for fixed and changing b 14|. However, the mean 
lifetime t(TV) grows exponentially with system size TV 
(Arrhenius law [bit ]) in the pure system case, making it 
stable in the thermodynamic limit. Instead, in the TD- 
CP, as b(t) can be adverse (b(t) < b CyPure ) for arbitrarily 
long time intervals, r(TV) is expected to be significantly 
reduced. But, does r(TV) still diverge for TV —> oo? i.e. 
does a truly stable active phase exist? 

Let us first report on numerical simulations of the TD- 
CP performed in dimensions d = 1, d = 2, and d — >• oo 
(for which we consider a fully connected network (FCN)). 
b(t) is independently extracted at each Monte Carlo step 
from a homogeneous distribution b € [bo — a, bo + er] . 
We performed both, homogeneous initial density experi- 
ments (with all sites initially active) and spreading ones 
(starting from a single active seed) [13| . In the first set of 
experiments we measured the average value, over many 
realizations, of the density of activity p(t) as a function 
of time and the mean lifetime r(TV) vs system size TV. 
In the second set, we measured standard quantities such 
as the survival probability as a function of time, P s (t) 
(l3j| . Searching for power-laws of the form p(t) ~ t~ e 
and P 8 (t) ~ t~ s , we determined the critical point lo- 
cation 6 ,c = 0.907(2), 0.656(1), and 0.500(1), and the 
exponent S rs 0.10(5), 0.126(2), and 0.5, for dimensions 
d = 1, d = 2, and the FCN respectively (in all cases 
rs S; a = 0.55 for d = 1 and a = 0.4 otherwise). For a 
fixed value of a, the shift &o,c — b CiPure is larger in d = 1 
than in d = 2, and vanishes in d — oo. Except for the 
mean-field value, and in agreement with previous findings 
11 1, these critical exponents are non universal, as they 
decrease upon increasing the noise amplitude a. Remark- 
ably, in d — 2 and d = oo, r(TV) scales at criticality as 
r - (In TV) 2 ' with z'(d = 2) = 5.18(5) (inset of Fig© and 
z'(d = oo) = 3.49(5) for a = 0.2. The values of z'(d) do 
not seem universal either, as they decrease with increas- 
ing a. Instead, in d = 1 we observe standard power-law 
scaling r ~ yv 1 - 55 ' 1 ). Furthermore, in d = 2 and d = oo 
(but, again, not in d = 1), we find a whole region within 
the active phase (b > &o,c) m which r(TV) grows generi- 
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FIG. 1: Schematic phase diagram for the pure contact process 
(CP) (solid line) the CP with quenched disorder (dashed line), 
and the CP with temporal disorder (dot-dashed line). For the 
second, a Griffiths phase appears within the absorbing region, 
while for the third a Temporal Griffiths phase appears (for 
d > 1) within the active region. The actual locations of the 
critical points may depend on noise intensity and dimension. 
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FIG. 2: Main: Log-log plot of the lifetime r as a function 
of system size TV for the TD-CP in d — 2 for various bo and 
a — 0.4. There is a finite region, bo £ [0.656, 0.675] with 
generic algebraic scaling of t(N) and continuously varying 
exponents. Inset: log-log plot of r(TV) vs ln(TV); from the fit 
at criticality (dashed line) we estimate r ~ (In TV) 5 ' 18 ' 5 ' . 



cally as a power-law with continuously varying exponent 
C, t(TV) — AK, with C -> as b -> b^ c (observe, in 
Fig 121 the slight downward curvature in the log r — log TV 
curves at criticality, reflecting the asymptotic logarith- 
mic behavior). Let us remark that obtaining data for 
larger sizes and deeper into the active phase, where the 
surviving times are huge, becomes excessively expensive. 
Hence, estimating with accuracy the upper limit of the 
algebraic scaling region is prohibitive. We have also mea- 
sured t(TV) in the absorbing state; as in the pure model 
case, it scales as r(TV) ~ ln(TV) in all dimensions. In 
summary, while the behavior of t(TV) in d = 1 is similar 
to that of pure systems, in d = 2 and d — > oo, we found 
(i) logarithmic scaling at criticality and (ii) an extended 
region with algebraic scaling. 

Let us now present analytical calculations for the high- 
dimensional limit (FCN). Given that, at every single 
step, the change on the global density p is ±1/TV, one 
can map its dynamics into a random walk in the in- 
terval [0, 1], with jumps ±1/TV occurring with probabili- 
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ties b(t) p(l — p) and [1 — b(t)] p respectively. The Mas- 
ter Equation for this process is easily written [Tij . and 
by performing a \/N expansion one readily obtains a 
Fokkcr-Planck equation whose (Ito) Langevin equivalent 
is (up to leading order), 



p{t) =ap- bp 2 + ay/pr](t) + (rp£{t), 



(1) 



with a — 2b — l + a 2 /2, b — bo, a — 1/VJV, and the noise 
£(t) = 2(b(t) — bo) /a. Observe the presence of both, a 
demographic noise, proportional to ^/p which vanishes in 
the N — > oo limit and an external or environmental noise, 
linear in p 0, [l5| . Generalizing to any spatial dimension 
one can show that the corresponding Langevin equation 
is just that of the directed percolation universality class 
[1 31 ] with a fluctuating linear-term parameter: 

p=(a + a£(t))p - bp 2 + V 2 p(x, t) + 7 \/^(x, t), (2) 

where 7 is a constant. For the sake of generality, we have 
numerically integrated Eq. (E|_in d — 1, d = 2 and in a 
FCN, by using the scheme in |l6| . We reproduced all the 
findings above proving that our conclusions are robust, 
and apply to any model in the directed percolation class 
not just the TD-CP. 

In the case of uncorrelated white noise (as used in the 
numerics above) the quasi-stationary solution of Eq. (fTJ) 



is P(p) 



P 



(1 + a 2 p/a 2 



exp(-2bp/a 2 ) 



which in the large N limit can be approximated 



by P(p) 



„2(a/cr 2 -l) p -2bp/a 2 



exhibiting an a/ a 



dependent singularity at p = 0. In this limit, the sin- 
gularity is not integrable for a < a c = a 2 / 2 (b < 
bo, c (FCN) = 0.5 for all a), and the only solution is a 
delta distribution at p = 0, i.e. the system is absorb- 
ing. Instead, for finite N, there is a 1/p singularity for 
any value of a and the only possible steady state is the 
absorbing one (as occurs for any finite system with de- 
mographic noise [Hj])- Defining z = lnp, Eq. (fT]) be- 
comes z = a — 6exp(z) + cr£(t), with a = a — er 2 /2, 
which describes a random walker trapped in a poten- 
tial V(z) = —az + bexp(z). It exhibits the three fol- 
lowing regimes: (i) Active phase (a > 0): the time 
required for the active state to fluctuate and reach the 
vicinity of the absorbing state (which is approximately 
p = l/N, i.e. z = — ln(iV)) and, eventually, die, is ex- 
ponential in the height of the potential [14|, t(N) ~ 
exp [V(-\nN)/(a 2 /2)] ~ exp(2a In N/a 2 ) ~ TV 25 ^ 2 . 
This is, t(N) exhibits generic algebraic scaling with con- 
tinuously varying exponents Hence, the active phase 
is truly stable when N — > 00. (ii) Critical point 
(a = 0): For sufficiently small values of z we have a 
free random walk (no potential barrier to be overcome) 
which covers a typical distance yfr in time r; equating 
this distance to z = In N, the time to die scales logarith- 
mically, t ~ [ln(7V)] 2 . (iii) Absorbing phase (5 < 0): 
z decays linearly in time and, hence, the time needed 



to reach z = — In AT scales as t ~ \nN. These predic- 
tions are in excellent agreement with the corresponding 
numerical results for the FCN. 

Result (i) can be recovered by using the path-integral 
representation of Eq. (fT]) 13j. The most probable path 
to the absorbing state can be easily calculated in semi- 
classical approximation, t is simply the inverse of the 
probability weight associated with such a path. By using 
this formalism, Kamenev et al. [l2j have recently in- 
vestigated in an interesting work the effect of correlated 
temporal disorder on a one-variable birth-death process. 
They conclude that, in the case of interest here (short- 
time correlated noise), t(N) grows exponentially with N 
for weak noise amplitudes and algebraically in the strong 
external-noise limit. Our result differs slightly from this: 
given that the demographic noise amplitude in Eq. (UJ) 
vanishes in the large N- limit, the strong- noise limit does 
not need to be invoked to obtain algebraic scaling. 

In order to extend our conclusions to finite dimensions 
and to make a parallelism between the reported broad 
regions of generic algebraic scaling —that we call tempo- 
ral Griffiths phases (TGPs) —and standard GPs, let us 
sketch the main properties of GPs for the contact 
cess equipped with quenched disorder, i.e. b — > 6(x) 
In the quenched CP, rare regions with b(x) > b CtPure 
and arbitrary size s appear with probability exp(— as), 
where a is a disorder-dependent constant. Such re- 
gions are locally active and, hence, activity survives on 
them until a coherent fluctuation kills it. This occurs 
at a characteristic time t c (s) ~ exp(/3s) where f3 is a 
constant, as given by the Arrhenius law [lij ]. Hence, 
the time-decay of the survival probability of a homo- 
geneous initial condition is given by the convolution 
P s (t) oc J ds exp(— as) exp(— t/t c (s)), and the leading 
contribution in saddle point approximation comes from 
size s*(t) = (l/p)\n((3t/a), implying P s (t) oc tr -!^ 
(right at the critical point, the exponent vanishes, and 
there is "activated scaling" , characterized by a logarith- 
mic decay P s (t) ~ (lnt) _e Q). Similar expressions ap- 
ply to the time decay of other quantities such as the ac- 
tivity density, as well as to many different systems with 
quenched disorder Q. 

Thus, some analogies between GPs and TGPs are: i) 
In GPs disorder is "quenched in space"; in TGPs it is 
"quenched in time", ii) In GPs rare (locally active) 
regions exist even if the overall state is absorbing; in 
TGPs rare (temporarily absorbing) time-intervals exist 
even if the overall state is active; i.e. the roles of ac- 
tive/absorbing phases are exchanged, iii) In GPs the 
probability for a (rare) active region of size s to occur is 
cxp(— as); in TGPs (rare) time intervals of length T are 
absorbing with probability exp(— aT); hence, the typical 
time to observe them is r ~ exp(aT). iv) In GPs, as we 
just argued, the leading contribution of the decay at time 
t comes from a rare region of size s* ~ \xt(tfi/a)/fi; this 
combined with (iii) leads to a generic power-law decay 
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FIG. 3: Log-log plot of the susceptibility, S, as a function of 
the field h, obtained by integrating Eq. |T} (main plot) and 
Eq. ([2]) for d = 2 (inset), for different values of a. 



in time, t a ^ . In TGPs, the time required to reach the 
absorbing state in an absorbing time-interval, is given by 
exp(-/?f) - 1/TV, or t* ~ ln(TV)//3. Equating t* with T 
in (Hi), one obtains a generic algebraic decay in system 
size; r(TV) ~ N a ^ . In conclusion, TGPs are analogous 
to GPs by exchanging the roles of space and time. These 
heuristic arguments seem to be valid even for finite di- 
mensional systems (down to d = 2). The reason why a 
TGP phase is not observed in d — 1 is not completely 
clear to us. Presently we are developing a semiclassical 
approximation, analogous to that in [l2[ , but for the spa- 
tialy explicit Eq. @ (see [3), m order to have a more 
precise understanding of low-dimensional cases. 

To further delve into the GP/TGP analogy, and given 
that GPs exhibit generic divergences of the susceptibility 
@j 1 we have measured numerically the susceptibility, 
defined as S = dp B t(h) / dh\h^o where p st is the average 
value of the activity in the steady state after introducing 
an external field h coupled to the system's dynamics, for 
Eq. (HJ) and for Eq. © in d = 2. Figfjshows that S mea- 
sured from Eq. ([1]) (which perfectly agrees with that in 
Monte Carlo simulations of the TD-CP on FCN) diverges 
all along the TGP. Fig[3]-inset shows generic divergences 
also for d = 2; however, given that very small values of h 
cannot be reached in this case, it is difficult to elucidate 
numerically whether such a divergence is for real in the 
h — > limit or it just a transient effect. 

To analytically understand these findings, we take 
Eq. (HJ) in the TV —> 00 limit and include an external field 
h. A new term exp [— 2/i/(per 2 )] , exhibiting an essential 
singularity at p = 0, appears in the stationary solution 
(see above). It is a matter of algebra to verify analytically 
that dp s t(h)/dh diverges algebraically as fo- 1 + 2 l a - el/< 7 ; 
in an extended interval a G [0,<7 2 ] around the critical 
point [15]. If, in this calculation, we replace a by 7 
which does not vanish when TV — > 00 and mimics what 
happens in finite dimensions (see Eq.(|2|)). the parameter- 
dependent singularities are replaced by the usual p _1 ab- 
sorbing state singularity. This suggests that the generic 
divergence of the susceptibility is a transient effect in the 
presence of non- vanishing demographic noise. In that 



case, the strong external-noise limit needs to be taken for 
the generic divergence to survive. Going beyond mean- 
field, it can be proved by using simple field-theoretical 
arguments (similarly to [15() that Eq. ([2]) with a = ex- 
hibits generic divergences of the susceptibility in a broad 
interval even in finite spatial dimensions. 

In summary, systems with absorbing states and fluc- 
tuating external conditions exhibit a region in the ac- 
tive phase —the "temporal Griffiths phase"— such that 
the mean lifetime scales generically as a power-law (with 
continuously varying exponents) of system size and log- 
arithmically at criticality. This occurs not only in mean 
field [HEH but also in extended systems as long as d > 2. 
TGPs have deep analogies with standard GPs, but the 
roles of space and time are reversed: in GPs (TGPs) spa- 
tial (temporal) disorder leads to generic algebraic scaling 
as a function of time (size). Moreover, as GPs, TGPs 
exhibit (at least in the strong noise limit) generic diver- 
gences of magnitudes such as the susceptibility and of 
stationary distribution functions. 

TGPs could be measured in the experimental realiza- 
tions of the directed percolation class with liquid-crystals 
[l9| by introducing externally changing fields, and could 
appear in many other systems such as in bistable Ising- 
like models with randomly changing conditions. We hope 
this work will stimulate new research along these lines. 
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